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ABSTRACT 


A general description of linear, switched-parameter 
systems is developed in terms of the mathematics of state 
variables. The restrictions imposed by the grouping of 
state variables and by the "core" states of the system upon 
the coefficient matrices are considered. Two procedures 
for determining the element values of those matrices are 
sae) Lllustrated with an example. A general expression 
for a Cost function to be used to measure system quality 
is developed and illustrated with two examples. Extensive 
recommendations for future work are made. Several examples 
Gos & ET Ce re eae Nie oeane 
for the improvement of the step responses of some switched- 


parameter systems are presented. 
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i. IN RRODUCT ION 


The first attempt to describe the behavior of a physical 
(real-world) system is usually based on a linear model of 
the system because such models are mathematically tractable 
and have been extensively studied. One purpose of this 
thesis is to develop a linear model using state variable 
theory for some real-world systems which may be modeled as 
imtcar systems with switched parameters. tTdhe parameters To 
mmewmruehned are the eoefiicients of the variables in the 
amavwe variable representation of the system and the values 
the parameters take reflect the actions of various switches 
iteeeae real-worla syStem. Altnougi thie meek Vive ie 
developed may be applied to some kinds of electronic cir- 
cumecs, this thesis deals primarily with linear, switched- 
parameter models of control systems. 

One example of the kind of real-world systems for which 
a model was sought is the arm positioning servomechanism in 
amagnetic disk-file. This servomechanism has two tasks. 
The first is to move the arm so that the read/write heads 
on the arm move from one "track" to some other in the 
shortest time possible. A "track" is defined as being one 
See ne marrow Concentric rings on the disk's surface which 
is used to store information. The second task of the 
servomechanism is to maintain the read/write heads (hence, 


the arm) in a nearly constant position with respect to the 








track of interest. The arm is positioned with respect to a 
reference set of pre-recorded tracks which reside on a 
special disk. 

During the transfer of the arm from one track to another, 
the servomechanism exhibits two modes of operation. While 
in the first mode, the system behaves as a "Bang-Bang" con- 
trol system, undergoing first a period of maximum acceler- 
ation, followed by a period of maximum deceleration. During 
the time that this mode is in effect, the position of the 
eee aS reckoned by counting tracks as they are passed. At 
the end of the deceleration phase of the first mode, the 
eaneis POSiItioned in the near vicinity of the new track. 
AG that time the servomechanism is switched to the second 
moeae wherein it benaves as a iilitar, nign-foin, Genrensesmeg. 
Geeesed—loop position control system. This second mode is 
Meeessary in order to obtain the required tracking accuracy. 
The transition between modes is accomplished by switching 
into the system a linear error detector and a compensator, 
and by disabling the "Bang-Bang" circuitry. 

The behavior of a system such as that just described 
may be modeled during each mode of operation by considering 
the “Bang-Bang" part end Une linear posi li Cmmeonmme! Servo 
part separately. However, the behavior of the system during 
Mme transition from one mode to the next is not describable 
by either of the two kinds of system alone. This thesis 


offers a solution to the problem by developing a general 


linear, switched-parameter model which is applicable to 








these kinds of systems and which allows one to deal much 
more directly with the transition between modes. The model 
developed uses state variable concepts and is reached as 
the limiting case of a general, ave order, linear, time- 
marrying system. The model is shown to be applit@ablesyo 2 
f@eooricted class of nonlinear systems. 

Previous work dealing with switched parameters in linear 
control systems has dealt primarily with the use of discon- 
tinuous naan iRaeron in various kinds of servomechanisms 
Odi. A different and particularly interesting application 
of switched-parameter techniques is that reported by Flugge- 
imteoez and Taylor [5]. The nonlinear system that they 
@escribe is a linear control system with multiple feedback 
Meems and a relay switching mechanism (oY Utveriiining water 
one of the feedback paths is to be used at a particular Time. 
The work presented in this thesis is believed to be the first 
attempt to develop a model applicable to the above kinds of 
systems. 

The second purpose of this thesis was to develop a 
flexible and general means of measuring the "quality" of 
a particular switched-parameter system. As an example of 
Weeat is’meant by the "quality" of such a system (and as 
another example of what is meant by "Switched-parameter 
system") consider the switched-parameter version of a velo- 
city control servomechanism which is to be used to bring a 
aioe tron a StOpped Condition Up vO. a =cervain Specit red 


movatvional velocity. The shaft is to be brought up vo 








speed by following a pre-planned acceleration curve until 
the shaft velocity is within the leck—-in vance sei soa 
phase-lock system which is to maintain the shaft at the 
@eerred velocity to within close tolerances. » ele teacuLs 

to the phase-lock system are a reference frequency and the 
output of a rotational-velocity-to-frequency transducer 
which is coupled to the shaft. At the instant when the 
rotational velocity of the shaft reaches a certain value 
faeein the eee range, the system is to be switched from 
meeacceleration mode to the velocity control servo mede. 
In such a system, one is concerned with how closely the 
programmed acceleration curve is followed, how much varia- 
Tion in the velocity exists while the system is in the 
Pemocity control servo mode, WriumwwewcAver vee. G8re ts 
Bme switching transients which may occur when the mode is 
Saeited. The "quality" of the system is some measure of 
how closely the system response approaches the design goals 
maeeach case. 

In order to have a single measure of the system's quality, 
the concept of a Cost function was taken from Optimal Con- 
trol theory and applied to these systems. Besides providing 
one number as an indication of the quality of a switched- 
parameter system, the use of a Cost function permits one 1G 
Meiize function minimization techniques in the design of 
the system. The use of such techniques was desired in order 
to “optimize" the performance of such systems. For instance, 


ose function minimization could be applied to the velocryy 


10 








servo just described in order to minimize the transient due 
to the switching from the acceleration mode to the velocity 
servo mode. Such a technique might also be applied to 
determine the optimum charge on one of the capacitors in 
the compensator which is switched into the arm positioning 
servomechanism in the disk-file in order to minimize the 
evavching transient in that system. 

The usefulness of the Cost function minimization technique 
is demonstrated fa the Appendix where the step function 
responses of several simple switched-parameter systems are 
"optimized" by determining an "optimum" switching time which 
[oon Zes an applicable Cost function. 

Preceding this demonstration is the development of the 


her 
3 Cm? oS} 


fi@e@el in section ll wnich 18S 1loLi@Weu, INT Secesen It 
the re deration of some constraints the physical world 
imposes upon the model and then by the presentation of two 
methods of determining the model from the real-world system. 
The two methods are illustrated by application to two 
different examples of the type of switched-parameter systems 
which have acceleration modes followed by position control 
modes. The applicable state equations have been solved on 
Pear eital computer to illustrate the behavior of the systems. 
fime Cost function 1s developed in secticoneiy and allus— 
trated with examples based on Cost versus switching time 


ferme the two examples presented in Section Tif. Section V 


is a presentation of a summary and some conclusions which 


uae 








are reached based upon the work done and the examples used 
to illustrate it. Section VI presents several suggestions 


ema recommendations of Topless fer itr geee ela 
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Ti. . MODEL CONSTEUC Hien 


the starting point for the constructionweone smote m a 
cal model to describe the linear, switched-parameter system 
was the general, linear, time-varying (LTV) system described 


as) tOllows: 


BaGo ) = ACE Coe (i ees 
N N N N (1) 
va) = C(t) xCr ot Geis 
~ ~ AY, ~ ~ 
Meee initial conditions, ie x(0), lee wakenea ies Om ibie- yas 


where 


x(t) is the n-vector of the system states, 
x(t) is’ thew time deri ya us7- mo. x(t) 

u(t) is the m-vector of control inputs, 

y(t) is the r-vector of output variables, 
A(t) ig an Nox nN Maurice ecOctn ne lemres 

B(t ) is an nm xX m Matrix Of cocifictemec. 
C(t) is ann x Y Matrix of (coetfiicients wane 


D(t) 1S a Ff XM Matweimomn coc tm) Cromer 


The determination of the As B, C , and D matrices is a 
straightforward procedure which is explained in many texts 
eae which will not be considered here. In fact, the deter= 
Manation of the coefficient matrices from the physical 
system will be less of a problem than may be indicated by 
the time-varying nature of those matrices in (1) due to the 


restrictions which are soon to be placed upon them. 


As 








Ilé is well known [6] that the time variation ef the os ace 
vector which is the general solution to the LTV system in 


GOe is given) by 


: 
x(t) = R(t)x, + X(t) Lo X(8) "BCs als ds (2) 


where X(t) iSMaAne Nexen matri< Satis tino ee following 
vector differential equation and where ae represents 
mifewmatrix inverse of X(t). X(t) is determined by solving 


the system of differential equations 
X(t) = A(t )X(t) 


Seer the interval te[0,t], with initial conditions X(0) = Zs 
where i is che n xX Heldenticy mauris 

Let the control interval [0,T] be divided into N sub- 
intervals of arbitrary length with the boundaries of the 


Subintervals denoted by the sequence of times 


where aN Ape | eerie > Ute for 0" < as he 
i - “i- = = 


ue 
Now let it be required of A(t), B(t), C(t), and D(t ) 


O 


that they all be matrices of constants for a Ss Se One 


fener o << min(t. -~ t.). i.e. the coefficient matrices are 
7 itl a 


constants except for: a small interval about the subinterval 


boundaries. If ts + 6 is denoted as oe and Us ee O sans 


denoted as eee then the set of differential equations 


describing the state of the system in the subinterval 


15 = : 
ve[t, 5 C547] may be written as 


Ge) ee ea aCe) BY baa ue) (3) 
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ia: we _ + + 
Weel Gnitial conditions cas x(t, ). The matrix A(t, ,t,,4) 
represents the constant value of the A matrix "we ener Paieuase= 
ular subinterval, with a similar interpretation apolicapie 


+ = 
ie B(t.t.,4)-. In the next subinterval, the A matrix value 


at: = 
c ), with a similar repre- 


would be represented by A(t. 45 $40 


Bencvation for the B Wiese 

Since the A matrix in any sarticular subinvery aia 
not, in general, have the same element values as the A 
matrix in any other Subinterval, and since the systems con- 
sidered are LTV, it is necessary that the A Mav Fix an Gene 
transition region between two subintervals be LTV. If the 
A and B matrices in the transition region are denoted hy 
A'(t) and B'(t) respectively, then the counterpart of (3) 


mememe Transition region may be WrivteNn as 
x(t) = A'(t)x(t) + B'(t)u(t) iy 
AV, AV) AV) %N Av, 


with initial conditions x_ = x(t~) and valid for te[t,,t.]. 
"NO av, al 1 1 


The form of the solution to this set of differential equa- 


tions is given by (2). The value of the states at t = t; 
is of particular interest since om) serves as the initial 
condition vector for the differential equations (3) which 


are applicable in the following subinterval, ee eee 

Then from (2) the expression for x(t) ) may be written as 
: 

mit, ) = Xte. )x(t;) 4 aC f_ ¥(s)7 7B! (s)u(s)ds 

= (5) 


where, as before, X(t) satisfies 


X(t) = At(t)X(t) 
™~ %™ % 


ite 








es m= + e ° ° e e e ar 
over the interval [t. ,t,] with initial condleions X(t, ) = 
~ ~\ 


mae the end point values of A'(t) constrained such that 


t a —_ 
A'(t,) = Alt, _1.¢,) 


> 
= 
Gir 
Ww 
I 
> 
= 
ct 


Poa 2 e 


Similar requirements may be imposed upon the Br Me Glee 
Since the interval [c,coe iS a transition perrod@ eur ing] 
which the system coefficient matrices A and B alter = varies. 
it is interesting to determine how the value of the state 
vector FG) varies as the transition period is made smaller 
and smaller. In particular, what is Lim x(t3)? From (5) 
mers seen that the interval over which the solution of the 
X(t) equation is required is just eye. This interval 


Gees to zero length as 6 > 0. Hence 


: re Es ae 
cae OG Oe 
Consequently, 


es 


a et iE - 
meet) = lim xe! )xe,) exec) ee eee 
S (ie IES py de he ee Qs) le ow ~ \ 
oo 
+ ; 
Using the limiting values of os and te and X(t.) > gives 


ee 
= lt — 
Tea Cero) ele aie a (a uence 
Sore Y 1 uy tL AY ts av AY av 


hee (6) 


emace the value of the integral taken over zero interval is 
Memo for a finite integrand. The integrand will be jfinice 


so long as u(t, ) is finite since Bt (t) is assumed to be 
~ 
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fmnite and since nde) a i as t; = a Hence saris Sr ae 
transitions between constant coefficient matrices are assumed 
to be no-elapsed-time switching of the coefficient values 
(which is accomplished by no-elapsed-time switching of 
system parameters) the values of the state variables will 
not be affected by the discontinuous change of the coeffi- 
ereme Matrices. 

pe far just the particular problem™er the ™transitaen 
between agilaeew, smiicieeciral within the overall: cenpcoum 
interval has been considered. However, the limit expressed 
mmo) now permits the entire problem to be treated in a 
straight-forward fashion. Because the differential equations 
(3) for each subinterval are known, the solution of (1) 
mene che piecewise constant iitweraerens on 4% 5, 
may be obtained by the sequential solution of N systems of 
n linear differential equations with constant coefficients. 


Thus, the differential equation 
X(t) = A(t. ,t,)x(t) + Blt, ,t,)ult) 


would be solved over the interval te[t ,t,] Wa Gly Giles setae ial 

conditions x = x(t_) = x(0). The value of the state vector 
GO 6 ~ 

at t = t,, x(t, ), would then be used as the initial condi- 


tion for the system of differential equations 
X(t) = A(t, ,t,)x(t) + B(ty,t,)u(t) 


which would be solved over the interval te[t,,to]- The 
final value of the state vector in this interval, x(t.) 


mowlda be used as the initial condition vector for theme. 


sl 


Subinterval. This process would then be repeated until the 

melution for the entire control inve vol siod been momm amor 
it has been assumed, to this point, that the vceoeiiaciens 

matrices are all altered at the same instant, at each of 

the subinterval boundaries. However, it is clear from (6) 

that the limit will be unchanged even if the B inal wee aS 

femaltered or if it is altered at some time either before 

ee atter the A matrix undergoes its transition. The case 

imor which the A and B Matrices undergo separate transitions 

meemconveniently viewed as a modification of the control 

out , u(t), caused by the variation of the B matrix, applied 

to a system for which the A matrix elements are switched at 

the times ty; or sees bya: 


oe es _ or are = - Aas oe SS ana ¥ fa ve 
Pme SYSVeMl may Ol COoUMS. HS Urcacea uC CUM nen 


t 


sequential, linear, constant systems by using some partic- 
ular expression for the time variation of the Al Tiel ae ak 
fering the transition period Ree subject to the con- 
straints on the endpoint values of the Al Mat rise Ment Vcied 
Meeviously. It is clear, however, that if the transition 
memerod 1s a small fraction of the shortest natural time 
@emscant of the system, then the instantaneously switched 
model of the transition will describe the behavior very. 
fame, lf the transition period is not small with respect 
to the shortest natural period of the system, then the LTV 
problem must be solved for the transition period in order 
mo, ObCain the initial conditions for the next subinterval 


mn which Che system has constant coefficients. The speed 
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of present day switching devices makes the no-elapsed-time 
mecansition a practical assumption for many systems. is.e 
mesultsS in the material that follows it is assumed Thaz 
meee Cranscata Onan. alison tne cect tiecrenitwme vice ce ee 
caused by an ideal (instantaneous) switching action. — 

For completeness sake, it should be mentioned that the 
Semcinuity of the state variables at the transition instants 
fees not imply that the components of the output vector, y, 
are continuous. Indeed, this is obvious from (1), since 
for x(t ) and u(t) DetheconPaanows LuUnectlons of tame, amy 


eescontinuity in the g or D matrices will appear as a 


smecontinuity in y(t). Let 
wv 


AC, = C(t, ,t,,4) 
and 


Se ee te een ee 


be the difference matrices for the switching at ts, Weer = 
1 < i < N-1. Then the vector Ay. = y(t?) =sy(t.Jawill be 
mee vector discontinuity caused by the discontinuous changes 


mimeme C and D matrices. Thus, 
nw NN 


_ + + 
ae = Clb, ote I Bty) + Bley st, ,) urs) 
a SAME aA) aM oo 


Having already assumed that x(t) and u(t) are continuous 
mmarections of time and, thus, do not change value at the 


eae cching instant, Ay. becomes 
VV 


Ay. 


ae ae (8) 


Ike, 








If, however, u(t) is allowed to change discontinuously 
(tC) is still constrained to be convinuens) eareenemine came 
Us > with 

Au. = u(t?) - Gee, 
vi Nv tL ~~ a reed 
then, substituting for u(t!) in (7), and with x(t,) = x(t7r), 
Nn ak N ali yw 


gives 


+ —. 
Ns SOc) i OCS a NOR) 


Hence, the difference vector becomes 


fore = ACox( t,t ADemers i) tect at. Au.. 
et Ci xh i? Dyul i? DC 5 ge i+1? al SP 
Mme lessor osf@@on (9) ger <né-Gieconkimuity cil wat) isqoeuwscoey, 
nw 


general in the sense that it applies to any discontinuous 
change of the Cc or D matrices or the control vector u(t) 
independent of any changes in either the A or 5 Mat eargecn 
A change in these latter two matrices will have an indirect 
effect om) since their change can be expected to 
influence the behavior of the state vector, x(t). 

The preceding work is founded upon knowing the answers 
MemeewoO questions. Both of the questions exist at the inter- 
face between the analysis and the design of the types of 
systems under consideration. That is, one must have first 
solved the synthesis or design problem before the two 
Questions may be answered and before the analysis may 


proceed. 
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The first question arises because the preceding deve lop- 
ment assumed 4 priori knowledge of the switching times 
ti> U5, eas tye: Because the manner in which these times 
are selected can completely alter the character of the 
mrvem, the effects of the selection methods snmeuld be con- 
Sidered. If the N-l switching instants are selected in 
some fashion prior to the beginning of the control interval, 
then the analysis of the system conceived as a sequence 
emeconstant, linear systems would proceed as described 
above. The examples presented in Sections III and IV and 
Appendix A are all of this type. On the other hand, the 


Switching instants may not all be fixed prior to the begin- 


meee Of Che control interval, but may instead be chosen 


with their selection being based upon some criteria involv- 
migetne current and projected state of the system. The 
real-world systems described in the Introduction are all 
examples of systems which determine the appropriate switch- 
ing instant by comparing the system's current state with a 
predetermined state which is fixed by the system designer. 
Mnese systems are intrinsically nonlinear. The nonlinearity 
occurs since the functional dependence of the switching 


instants must be displayed as 


bt, = ty (x(t), ult) ,t). 


This functional dependence of the switching instants requires 


that the coefficient matrices be considered, in general, 


a 








mov only functions of time, Dut Tuncttonewor the wore ome 
eentrol vectors as well. 
The functional dependence of the A and B matrices may 
QW 


be displayed as: 


(Cr ewer). ) ; 
_ om excep 
(x(t) ,u(t),t) 


CW e8LS 
T 
fe Coe 


with similar expressions for C and D. That the system is 
mamerently nonlinear is evident since the character of the 
system's response will be, in general, dependent upon the 
imouLv. However, the matrices are still constrained to be 
Pepe cCewise constant over the various subintervals of the 
Pomatrol interval [0,T}]. Consequently, the general functional 
relationships disnlaved above describe only formally the 
mmaievional dependence of the very narrow class of systems 
wWaach are considered here. . 

ifhiemonce the aswitehang@sinstant -_ssdevermined .aume 
transitions of the coefficient matrices proceed without regard 
mem the state or control vectors, then the limiting argument 
used in obtaining (6) shows that the states are still con- 
Manuocous at the switching instant, and that the linear, 
piecewise constant model can still be used to give precise 
results for what is technically a nonlinear system. The 
applicability of the model to a variety of real systems is 
meus assured. 

The second question that requires an answer before there 
can be a completed mathematical model probably has no 


"right" answer. It requires that one have determined the 


ae 





best structure for the system which is to De employed. if 
the control interval is divided into two or more non-zero 
subintervals, then it is necessary to determine which 
sequence of coefficient matrices is the best. Alternatively, 
one may seek that sequence of linear, piecewise-constant—- 
memeamneter control systems which would best cosine the 
required task. This problem may be formalized as that of 
Pending, for instance, the A matrix Of “a particular scyve tem. 


ime. find the function, fi; such t hee 
To) 


A(t, oti) = f(x(t) u(t), t,a(t,_45¢,)) 


Meme <= 1 <= N-1, will cause the resultant system to have 
the desired characteristics. This problem remains to be 
Zeno,  68ie OCSL Vitae Mae Ivo em lCccec, St resent, «* “aw 
meme SsOrt of heuristic approach which is aided by the 
physical constraints that a particular application imposes 
will lead to a system which gives acceptable performance. 
From the preceding it has been seen that linear, 
Switched-parameter systems are easily handled by the mathe- 
mewics of the state variable theory. It is apparent, how- 
mer, that the set of coupled differential equations that 
describes one of these systems is probably more amenable 
CTO being solved on a digital computer than “6 being ATER’ 
in a mathematically closed form. The closed form expression 
is possible in principle but would become increasingly 
clumsy as the system order and the number of subintervals 


increased. 


es 








Another possibility is that hybriG@ commuter toca. 
might be used to obtain a solution to the state equations. 
im applications that do not require @ great cealwor 
accuracy in the solution, such an approach may prove useful; 
however, the insensitivity of the hybrid computer to small 
mie, Significant changes in a parameter value or the length 


of a control subinterval may not allow the use of these 


techniques. 
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Iti. APPLICATION TO PARTICULAR LINEAR Soo Cis 


ee LE ORE SRT CTE BSS RUCTURES Cr Shih. eo) 

The general switched-parameter system which was the 
concern of the previous section will seldom be encountered 
in its full generality. One reason for this is that physi- 
memesystems exhibit apart of their reality as restrictions 
on the forms that the coefficient matrices of the state 
variable representation are free to take. As an example of 
this sort of restriction may be cited that shown by the 
@rassical linear description of a servomotor. The Lean tet 
Mmmrrection iS shown with a block diagram of the servomotor in 


Pivure i. 


u(t) —----—_ y(t) 
Pipe a ServoncLor= bihecl Daler, 


The second order differential equation which describes the 


behavior of this servomotor is 
ite) se eyeGr)) = lepiee))e Gee 


where the dots indicate the time derivatives and a and k 
eme constants. 
ie puUby IIe) imMbowsbaverVvarlable. form aves —eOllmon 


practice to choose as the states of the system the function 
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w(t) and ats time derivatives. Um eo omllcr 


Cece 


and 


i GD 


be the system states. Then, rewriting (11) in state vari- 


able form: 


x, (t) BAG e) 


: (12) 
X5(t) = -ax,(t) + knew ) , 
which may be Dic Maen CiOnsinal cies) Orlnwas 
ee OF Lite. Ce) 0 
ads tte] tfuce)), (13) 
x, (t) QO -a x5 (t) k 


CNG 


role SeeaCemes 1p Ge) e 
nN" NY nN 


The point of this elementary exercise has been to show that 
the forms of the pean the B matrices are restricted by the 
particular device and by the quantities selected ae cue 
mmeve Variables. In (13) the elements of the first column 
of the A matrix will always be zero for this particular 
emorce of states to describe this particular device. Addi- 
tionally, the element of the B matrix that is presently 
Zero, Will always be zero under these conditions. Indeed, 
one may say that for an nee order system which has this 


servomotor in it and which uses y(t) and y(t) as the two 


states that describe the behavior of the servomotor (and 
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n-2 other states to describe the rest of the system), if 
x; (t) corresponds to y(t) and x(t) corresponds to y(t) 


(i # j), then 
x; (t) = x(t), (14) 


meee ene Only non—=zero Term ani the oe POV Ot Ste Sijcmme mas 

2 and z MSE eI Oa) Valk love ass =e ite Pes v raC mse les © tee are 
elements of the two coefficient matrices result because the 
fame omOovor 1S deseribed by a second order dil ferenvial 
equation which, when expressed as coupled first order dif- 
ferential equations will always result in (14) as the first 
order differential equation for x, (t) (= y(t)). Furthermore, 
the two states x, (t) and x(t), or their equivalents, will 
alwavs be present and always be paired as in (14) in any 


system in which the servomotor is used (still assuming the 


basic y(t) and y(t) choice of states). 


Pee lie CORE STATES 

The concept of grouping the system's state variables as 
occurred with x, (t) and x(t) above forms the basis for the 
TOllowing discussion of switched-parameter control systems. 
No matter what the system's order may be, there will always 
Seema 'core” of states which one is trying to control. For 
instance, ina shaft positioning servomechanism, one seeks to 
eontrol the states of the servomotor which drives the shaft 
eae those states (shaft position and shaft rotational velo- 
city) become the "core" states of the system. It may be 


that various compensation networks have been employed in 


any 








the system to improve its response. The state variables 
needed to describe the behavior of these networks, taken 
together with the "core" states, make up the n states 
necessary to describe the system. 

The essential point is that one seeks to control the 
values of only a few of the total number of the system's 
states. This does not mean that the remaining states are 
not to be controlled in some sense. It does mean that one 
is less concerned with their endpoint values than with the 
effect they may have on the variation of the "core" states. 

One is led naturally, then, to questions concerning the 
order of the variation of a particular state (or set of 


states) within a system. If the system can be represented 


Meroe Ol Titrocu Order Git tperventie: caustics cs 
em Ce) Ae ene Gt) B 
mages) ary EE a | ee ee [ue Gi 
ey 0 A x B 
16 eer x, (t) all 


ee) ; 
where, for ann order systen, ae 1S Gan oe Mane ese, AG 
“MV 


ieman r xX (n=-r) matrix, A, is an (n-r) x (n-r) matrix, and 


nl 
0 is a (n-r) x r null matrix, then, one may say that the 
states comprising the state vector x, (t) will exhibit varia- 
moms characteristic of an aa) order system since uney 
are totally decoupled from the r states represented by the 
x, (t) state vector. The states of the x, (t) state vector 


will exhibit the variations characteristic of an pen Srouer 


system, 


; 
xO {t) AQXo{t) : Rox, 6b) se sy eeae J) 
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where the last two terms represent the forcing function for 
the inhomogeneous system. 

The linear, switched-parameter systems which are repre- 
semcted in a particular time interval by (15) are 2 subset 
of the more general class of switched-parameter system 


described by 


S } t 

Te. Z BO, ofa as : Fo [uct )| alae 
: 1 | 

X,(t) Ay 1 Ard Li 0*? Bi 


maemeen is the same as (15) exeept that the (n-r) x r null 
matrix is generalized to the matrix Aj. The x, (t) state 
vector for this case (as for (15) when that equation is 
interpreted as describing a switched-parameter system) 
represents those of the system states which are considered 
the "core" states. 

The switched-parameter system to be considered will be 
described by a sequence of N linear, constant-parameter 
systems during N successive subintervals of the overall 
control interval which is taken as [0,T). During the i%® 
subinterval the system is described by (16) with the initial 


Pemditions, x. = x(t.), and for te[t. ,t With A! and 
wo a Jt ad 


a 


le, taken as non-null matrices and assuming that A is rank n 
Nw “aw 


sag! 


(non-degenerate system), the "core" state vector, x, (t) may 
have ae Onowelsio lelouwivikone, dig Wall 1 Jere elats lilac elo sie waa Oe letse 
ime partitioning of the A matrix in (16) serves to delimit 
the "core" states. The system's differential equations have 
foam so —Orgdercasthac Such @ partitioning is pessible. “Neu, 


assume that in the interval [t ]) the behavior of the 


4410442 


Zo 








th 


moerce states 1s To. be less et noamen order. Such behavior 


could be caused by forcing the Aj, matrix to become [O:A,°] 


in the Gee subinterval. This may require that the 
equations of the states not in the "core" be rearranged. 


The discussion will proceed as well if the As Matra. 1s 
Nw 


merced Co become the null matrix. Thus, the behavior of 
the "core" states will be altered from ae order to Goole 


order at t = t,,,. It is understood that the "forcing" 


of the At matrix is to be done by switching the values of 


tme system's parameters. The next linear system to follow 
that described by (16), then, will be the system (15) with 


) ame valkia fomete lt ioe). 


a] Diels sige 


tige initial conditions for (15) may be viewed as the initial 


metal conditions x. = x(t, 
NO rand 


Seg@encolons Loretire Tiere suvinvervel (xéc)) tremstormed wv 
fv 
The intervening sequence of linear systems to their present 


meme, x(t, 
nJ 


t+)" The effect of the preceding systems is, 


thus, apparent. It may be noted that if knowledge of the 

variations of x(t) is not required for further work, the 
{iV 

em@eare system of n differential equations need not be solved. 


One may solve for this special case just the equations 


X(t) oo i eee (17) 


with the initial conditions being the appropriate subset os 
@@e initial conditions for (15). 

Since the order of the behavior of the "core" state 
vector is (n-r), which is assumed to be the number of the 


meere states, the system will be restricted in the Changes 
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that may be made in it at the end of the current subinter- 
val. That is, the character of the behavior of the "core" 


states may remain of the same order, (n-r), or the order may 


ne order) by suitably altering the A, 


and the AS matrices to couple additional differential equa- 
nr 


_ 


be increased (up ton 


tions to the equations of the "core" states. If the order 
is to be increased, the values of those states in the x(t) 
vector which are to be coupled to the x,t) state vector 

at the next switching instant will become quite important 
Since their values at the endpoint of the current subinter- 
val will be used as the initial conditions for the set of 
Sempled equations which include the “core” state equations 
in the next subinterval. Hence, one will be solving (17) 
e@emented Oy soli @aMittenal stabe@, ~'(t), where tne 


" matrix 


coupling is accomplished by the equivalent of the A, 
mentioned earlier as the non-null part of the A Mae 25a: 


i.e. one must solve 


x} (t) = A,x,(t) + Auxs(t) + Bru(t) Cdlis)) 


where the initial conditions are given by X1(ts44) and 


t t 
xh ae P  toOretne ssUbanter ra | telt,,oteyol> and where Btu(t) 


represents the appropriate forcing function terms from (15). 


it] 


Ptetiay De WOvea sin passing Thaw thevinitial Conditions 
represented by xi(tsay)? may be controlled to some extent by 
Che system designer during the subinterval in which they are 
decoupled from the "core" states. This flexibility is 


allowed since the vector of states, x (t), in (15) may be 


Sil 








treated as an independent system so long as none of its 
member states are coupled to the "core" state vector, x, (t). 
The improvements in the response of a switched-parameter 
system that such flexibility may permit remain to be deter- 
mined. The values of xh(0) (for the case in which low order 
behavior of the "core" states is desired during the first 
subinterval) are a particularly important example of the 
fexability that is available. 

Once the linear model of the system has been placed in 
the form of (15), and the elements of the A matrix that are 
to change and the values that they are to take after each 
change have been established for each segment of the control 
maperval, then the solution of the system's differential 


e '_ 
ee ee ee | W TT CN WITT 
_ ww 


eamaticns May Preceec strakent 


C. DETERMINING THE STATE EQUATIONS 

mis obvious from the foregoing material that the 
Meseription of switched-parameter systems is most easily 
mecomplished using state variable concepts. Furthermore, 
miuetne state variables are restricted to those which have 
an immediate physical significance (e.g. capacitor voltages, 
mauccOor currents, shaft rotational velocities, and the 
like) the conceptual understanding of the switchings of the 
system's parameters and the resultant behavior of the 
system's states is enhanced since one can readily see how 
meh of one particular state is being coupled to some other. 
The order of the system and the values of many of the system 


parameters will have been fixed during the design of the 
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system. Usually the fixed parameters will be fixed because 
of the inherent physical characteristics of the system's 
components; although, some parameters may be fixed by 
"designer's choice." Those parameters that do not have 
fixed values during the entire control interval for one of 
ie above reasons are those that are to be switched in value 
imeem subinterval to subinterval. 

Two procedures present themselves as possible methods 
WONG translating the physical system with switched-parameters 
mmtvo the state equations of (16). The first procedure would 
have one determine the coefficient matrices CA, B, C, and D) 
assuming a maximum amount of coupling between the state 
Semations; that is, utilize individual multiplicative 
MaAViN=e & Wolwe CO ts L 
the switching functions are non-zero, and determine a general 
set of coefficient matrices. Then, by allowing the appro- 
Priate switching functions to take their zero values in The 
proper subintervals, one may read the element values of the 
€@eificient matrices which are appropos of that subinterval. 
This procedure requires that one manipulate the differential 


ne order system (assuming that the maximum 


equations for ann 
order of the -system is n), complicated by however many 
switching functions are required. for the case of N sub- 
intervals, there would be N-1l individual switching functions 
necessary. 


The second procedure asks that one proceed on a sub- 


interval-by-subinterval basis. Beginning in the first 
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subinterval, the state equations GbeunCse sso ea Geom i eetiee 
coupled to the "core" states are written down. Next, the 
equations for any other states which are not coupled to the 
"core" states but are coupled to a forcing function or 
other states not in the "core" (but are not isolated) are 
Peer Len down, leaving only the iseolaved Svaves Unaccounvea 
for. One then considers the system as it is to be after 
(ime first switching instant, identifies the isolated states 
which have been coupled to the "core" states, and adds rows 
ome. columns to thes appropriate system coefficient matrices 


Mai speigieis Tee claconois  — JMcre 


mamorder to take these 
esroments of the coefficient matrices that correspond to the 
"old" states are changed appropriately so as to take account 
Be viie COUpLITE Of tte @ivoupliuB Of uNMOSe sudtes. Ge CUR 
tinues in this fashion, adding rows and columns (adding 
differential equations) to the coefficient matrices until 
they have attained their maximum dimensions, after which 
time the coefficient matrices are determined simply by 
altering the appropriate element values for each subinterval. 
Care must be taken to ascertain the ial value of any 
"new" states when they are switched into the system from 
their previous isolation the first time. 

There are some pitfalls to be avoided in both of these 
meocedures. For the most part, the pitfalls are related to 
the identification of the states of the system. In switching 


parameters, it may occur that one obtains the equivalent of 


meer Culm which there ismasloop containing all capacivers, 
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erea node having only indwetrors conmectcd aren bam eu 
system is degenerate in that not all of the "states" repre- 
eemcred by thescapaiemtor voltages, OMbhe Sigducter currents, 
ere independent. This will require that the A Maids ois 
the system in that subinterval have rank one less (assuming 
only one degenerate state) than it would have if the system 
were not degenerate. The system of equations may still be 
solved but methods that require inversion of the A Menuet X 
must be used with care. The USC Gileeberativesdal tenrential 
equation solvers such as the Runge-Kutta and the Hamming's 
fee@iiied Predictor-Corrector methods will provide time 
meeponse Solutions to the state equations with a minimum 
Meeeconcern about the values of the various coefficient 
iaerices. 

Another type of pitfall may be encountered when using 
the first procedure mentioned above. The problem occurs 
when allowing a switching function to take a zero value 
meres an Clement of a coefficient matrix to become infinite. 
This is probably the result of dividing by a switching func- 
tion which was being treated as a symbolic, non-zero quan- 
Mey during some algebraic manipulations. The problem is 
best cured by guarding against such manipulations. 

One other point needs to be mentioned. The system, in 
each sub-interval, appears as one or more sets of coupled, 
iemear differential equations. In many cases the system 
Being modeled will be a classical, linear control system 


With a single input and a single output. Such systems are 
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well described using transfer Funct tens, pole -Zerenr lou. 
and the other well-known methods applicable to classical 

Ss svems. Unfortunately, these classical eanalysisegecr— 
nigues are not well suited to the description of switched 
parameter systems even in the single-input, single-output 
case since the identity of the individual states is easily 
feoev. ilhis is not so important during the period that the 
eascem 1S in one particular subinterval, but it is very 
important at the ends of the subinterval because the values 
Seeecone individualstates must be known in order to determine 
eee anitial conditions for the subsequent subinterval. For 
instance, if in one subinterval a system has a transfer 


function, T(s) which changes as shown: 


ol ae 
and 
T(s) = 1 for >t 
2 ates 
s (sta) 


time Dehavior of the system's output will depend on the 
Semucture of the second system. If the system in the sub- 
mirerval prior to t. is structured as three individual l/s 


terms, then the system which results from the switching of 


(2)(2) ~ (sta) (2](3} 


will have a different behavior than if the switching goes as 


E12) -BIEI es) 


i 


alt 
Ss sta 
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This will be~so because the initial condi Gicneonm uno) 
term will vary with its position in the system and thus cause 
the transient response term of the system output Co vary. 
Since one must look at the individual state variables in 
order to determine the structure of the system in two adja- 
cent subintervals, it is convenient to describe the system 
@@ey in-terms of the state variables. Of course, it is 
me | possible to consider the transfer function for a 
Derticular subinterval where that is possible and where it 
may lend additional insight into the behavior of the system 
in the subinterval. This will usually be advantageous when 
the effect of the parameter switching will be to lsrealiayee alg Ne 


to isolate an entire group of states such as those describ- 


A simple example will be worked for each procedure to 
illustrate its use. The two examples have been selected 
Pteemathematical simplicity and as illustrations of systems 
that show improved step response over purely linear systems. 
Both examples are third order systems with only one switch- 
ing of the parameters considered. Single-input, single- 
output systems were used with the output being one of the 
states in order to simplify matters even further and to 
remove the necessity for considering the C and D Weer ieces 
in the examples. Because of the diagrammatic simplicity, 
the systems are shown in a block diagram form in the time 


domain with each of the states represented by an integrator. 


Sti 











temo ,'l'] 


Figure 2. Block Diagram of a Third Order, Switched-Parameter 
Control System which has a Core of States Showing 
Third Order Behavior Before and Second Order 
Behavior After Switching. 


Figure 2 shows the system that is used to illustrate 
merireavion of the first procedure. “ne symboi Sy represents 
the switching function and takes on the two values indicated. 
The identity of the states, as shown, is arbitrary provided 
Only that they remain distinct. The first procedure requires 
mice all of the switching functions be considered non-zero. 
This is easily done by treating the functions (in this case, 
only S,) symbolically while writing down the n coupled 
Gifferential equations that describe the system. Thus for 


the system of Figure 2, one writes 


x, (t) = (Mase) 


X(t) X(t) 
X(t) = Ss, u(t) + (1-s,)x,(t) + (l-s,)x2(t) 


48) X3(t) - s,x,(t) 





or 


0 0 0 CS ) 
x(t) = 6 0 1 x(t) ce Ge) 
VSS) So elas] adem Sy 
mor Sy equal to zero in the first subinterval and equal to 


mae in the second, and with the switching taking place at 
t = Ci> the systems on either side of the switching point 


are readily seen. to be 


0 0 0 al 
Sele ee a0 0 Ui lseGt eect | Ome Ct.) eons @kep' 
Vv ~ ~ 
1 0 iL 0 te[0,t,] 
and 
0 0 0 0 
x(t) = 0 0 Lix(pyes | 0] ue for (20) 
0 =i =u il te[t, ,T) 


meere the element values of the & and B New rt Ceseiavem eed 
developed explicitly. For this simple case the procedure 
Womeks well. However, when the system complexity increases 
and there are more than several switching functions to con- 
Sider, then the expressions required for the individual 
elements of the various matrices oe Kec anly wmamiveleily 
Phe second procedure should then be used in preference to 
taas one. 

The two sets of equations (19) and (20) were solved on 
a Gigital computer by means of an iterative routine that 


uses Hamming's Modified Predictor-Corrector method for 


oy 








solving sets of linear differential equations. The method 
utilizes a fourth order Runge-Kutta procedure to generate 
four initial points, and then proceeds, using the preceding 
four points in an iterative procedure to generate an esti- 
mace Of the next point, and, by extension, the Solution Go 
the set of equations. For further information on these 
methods see [7, 8]. Figure 3 shows the time response of 
the states of the system for two different values of the 
switching time, ae The input, u(t) was a 10 volt step 
mumebLOn in both cases. The length of the control interval 
mere O Seconds. The initial conditions of all the states 
were taken to be zero. The switching time for Figure 3(a) 


eee = 1.45 seconds, and for Figure 3(b), oy = 1.349 


il 
seconds. 
Figure 4 shows the time response of the states of the 


pemema Order System Wath constant parameters, 


: 0 il 0 
x(t) = x(t) + foc ele 
ae -l] -4|~ 1 
memo 10 volt step input. Zero initial conditions and a 10 


second control interval apply also to this system. Figure 4 
is intended to serve as a reference for comparison with the 
time responses of the other systems. As may be seen from 
the comparison of Figures 3 and 4, the switched-parameter 
eyovem shows a marked improvement in step response. 

The comparison of these two different systems is justi- 
fied since the switched-parameter systems may be viewed as 
being the same second order system with means for modifying 


Hes iiivaal COndivions. — Thus, one is coOilparing the reste 
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Figure 3. Time Response of the Third Order, Switched- 
Parameter System of the First Example to a 10 
Volt Step Input. ca) switch new lame. ty es 
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Figure 4. Time Response of a Second Order System with 
Constant Parameters to a 10 Volt Step Inpur. 
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of the constant-parameter, second order system when started 
from zero initial conditions to the response of the same 
system during the time its initial conditions are being 
"pre-charged" by the unstable third order system (in the 
first subinterval) and after it is released as a second 
order system (in the second subinterval). . 

The choice of switching time, ty> for Figure 3(b) was 
decided by repeated trial-and-error solutions of the state 
equations. The value obtained, t, = 1.349 seconds, gave 
the best-looking response. An increase or decrease in this 
value by as much as 0.001 seconds produced a visually 
detectable increase in the overshoot or undershoot respec- 
tively. A similar sensitivity to the value of the switching 
Mememtitas veer noticed in @11 Of Wine j=xampies wiFreu. Fite 
effects of such sensitivity may preclude the use of hybrid 
eomputer techniques in the sclution of the state equations 
of such systems. 

The second example will be equally simple. Figure 5 
shows the third order system in the first subinterval. As 
before, the states have been individually identified and 


iepresented by integrator blocks. 





EO ste a 


eeeure 5. Block Diagram of the Third Order, Switched- 
Parameter System of the Second Example Prior to 
SVL ella lak @ yea 
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fhe system differential equations in this subintervalyare 


0 0 0 il 
(ey = 0" on eee) eae uct} , (21) 
lad ra) 
il 0 0 0 
te{0,t,]. 


It is desired that the system take the configuration shown 
maeeagure 6 in the subinterval te[t,,T]. Ehe weiam =e.c ant © 
be made can be systematically determined by examining the 
Seerterential equation of each state in turn, as they are 


Baven in (21). 


u(t) 





we oes 


Figure 6. Block Diagram of the Third Order Switched- 
Poaremeter yoysvem of the second pees ee 
Si bela Mees - 


moerae from Cop to bottom of the state vector, examining 


x, (t), X(t), and X(t) in turn, it is apparent that the 


new set of state equations becomes 


yd 
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0 0 0 0 
x(t) = 10 0 DP ssh) S| 6 fact) . for > 


Oo -1 -4 i. 


If an isolated fourth state had been encountered in the 
Beeond Subinterval, and if the coupling between the states 
x, (t), XK, (t ) and X(t) Was as ShOWno in (ec) sec e qeariae | Moe tla 


Sigeer System of this subinterval would be formed as 


, 0 
| 
x(t ) 0 0 NR Gp) 0 
———|= ae | = AF fut )| (23) 
rete) 0 -1 -4Yo Kole ) 1 
a a ee gb 
Hal : Ay al 


ee Uh m1 
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where the elements of Ay > A,> AS? and Bo would be determined 
from the coupling between x(t) and the three prior states. 
imesmould be well noted, that since this is the first time 
that X(t ) has been used in the system (else it would already 
be included in x(t)), that particular attention must be paid 
to the value of x,(t,), which is the initial condition of 
x, (t), Lometne mevmolloim tei yada 

Depending upon the method used to solve the set of 
differential equations (22) or (23), the chore of program- 
meee a digital computer for the task of solving those 
equations may be eased by always working with an oe Order 
system, even when there are isolated states that are being 


"ignored." It would then be necessary only to alter the 


oe 








element values in the coefficient matrices and not necessary 
to change their dimensions during the course of the problem 
solution. The isolated states are conveniently allowed for 
even if not individually identified, by adding the equiva- 

t 
1° 81° r 


some trouble encountered in determining the maximum system 


lent of Ao ok and Bo as null mat ries es i eine cma: 
order, the lower order system's equations may be written 
down and then adjoined to the appropriate null matrices, 
when the maximum system order is known. 

As for the first example system, this system, as des- 
eabed by 21) and (22) was solved on the digital computer. 
The trajectories of the states are shown in Figures 7(a) 
and (b) for a control interval of 10 seconds and a 10 volt 


bape Ses 
{ 


9 = e 
= = ~~ Ce lone hd Taare 
ty a cwi tons 


So eon hinc Tages 


eeoo LUNCtLON inpuv. Freure 


t. = 1.349 seconds, and 7(a) is for t, = 1.699 seconds. 


il il 
Figure 7(a), when compared with Figure 3(a), shows that 
the response of this system has a slightly longer rise time 
than the first system for about the same amount of over- 


shoot. The response of this system is also markedly better 


than that of the plain second order system. 


46 








16, 





LiJ 
QO 
=) 
z= 10 
=} 
ae 
> 
ey < 
5 
X3 
Ne ED 
O ti2 4 6 8 10 
TIME IN SECONDS 
af 
ae | - 
AY 
WJ 
Q 
= 
= 10 
—) 
5 
foe < 
5 
Pieces lian inaaalye ayeeNe 
O Lia? 4 6 8 10. 


TIME SIN SECONDS 


mreure {. Time Response of the Third Order, Switched- 
Parameter System of the Second Example to a 10 
Volt Steps bapuc.me (a) SWitehune. lames ty = 1.699 


Seconds. —) (bd) ty = 1,340 seconadc. 


af 








iV. MEASURES OP gS is EU el 


The various methods available for taking the measure of 
a system (or its model) have the common goal of providing 
the analyst or designer with a quantitative estimate of how 
well the system may be expected to accomplish its designed 
tasks. Linear system models are commonly tested by applying 
a known input (usually a step function or a sinusoid) end 
Comparing the output with the desired output. In the case 
of the step function input, one might examine the rise time, 
o@gemsh)ool, and settling time of the output response. As 
may be seen from Figures 3 and 7, there seems to be justi- 
fication for using tnose measures of system Gualivy to 
grade some switched-parameter systems. 

Unfortunately, the usefulness of the classical system 
response measures is somewhat degraded since it is usually 
not possible to make the direct correlation between the 
classical response measure and a system parameter that is 
possible with the classical linear system. The Bene eeien 
results from what may be considered an inherent ambiguity 
in linear, switched-parameter systems. The amOiguity arises 
Meese there may be more than one sequence of linear, 
constant-parameter systems that would have the same output 
moreade Piven inpuc.  PWhis is another manifestation of the 
intrinsic nonlinearity of switched-parameter systems. Still, 


for those switched-parameter systems that are designed to 
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replace a classical linear system, the classical response 
measures may well provide a useful means of grading the 
systems. 

There are cases, however, for which the classical 
Meeoures of system performance are not adequave- As one 
example may be cited the Linear Regulator problem of Optimal 
Control theory. A Cost function approach to the measure- 
ment of the quality of linear, switched-parameter systems 
has several features to recommend it. The paseiii thes of 
having several subintervals in which such systems may have 
Breatly varied behavior and in each of which the designer 
may. have separate goals, strongly recommends a flexible 
Pees ure Oe com Quad ty. ny addiGlon 4 Oster unel tom 
Peeve use in PROWL amiieasure of sewterai variatlies 
in the system at the same time. For instance, one might 
construct a Cost function for a system whose value would 
be the combined measure of the system's settling time and 
Che integral of the error between the output and the 
input squared (assuming a system for which such things 
make sense.) Another reason for desiring to use Cost 
functions with switched-parameter systems is that such an 
approach lends itself to some sort of system optimization 
mecomolished through minimization of the Cost function. 
Several examples of Cost function minimization are presented 
Wms pendix A as a demonstration of the utility of the 
method. The examples are limited to the determination of 
the optimum switching time for some simple systems but they 


illustrate a method that might profitably be employed. 
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Against this potential flexibility must be balanced the 
loss of detail that is inherent in the use of any sort of 
Cost function. The loss of detail arises, of course, 
Since the Cost function is the weighted sum of several 
meberent quantities and there is always the possibieicy 
that an increase in one quantity was offset by a decrease 
maeenother. Ambiguity in the value of the Cost function 
mmpossible when the changes in the different quantities 
offset each other so that the function value remains 
unchanged. 

A general Cost function for switched-parameter systems 


may be written as 


Jm = QT, + RT, (24) 
Pier eC 

Q = [oe eee Onl and Q5 Pe Ooi ie lie 

ee 10) 


are weighting factors for the Cost vector, J.; and the 
taba of the control interval, oe respectively. The 


vector 


ah 
[J, J nN! 


Jy # Wy sTgs-5d 


the length of the control interval is defined as 


N 


oo, nee gir 
ie 1. 


for the switched-parameter system given by 
Ge) aoe sr) b siGa)wilGy)) 
nw ww Ww lad Ww 


mad Valid for t in the interval [0,7 J. 
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In the i°” subinterval, 


Wee 
Jy = wltg yy) Hywlta yy) #2 7" Uys) "Myy(s) + us) "Nyuls) Is 
a 
(25) 
ea ee PS sige) 


where He M., and N. are 1 xX Nn, n xX Ny and Mx en postEaye 
IY ws 


semi-definite diagonal matrices of weighting factors, 


igemyoectively. The vector, r,(t), is the desired trajectory 


iomextt) in the subinterval te[t,,t.,4], and x(t, .4) is 
the desired vélue of x(t} et t = tay z+ An equation analgous 


meme>) 28 applicable to each of the N subintervals. 

Mien of the weighting factors are to be non=Zzero, and 
what the values are to be is left to the system anaiyst 
MemeeaecciSion since the weighting factors should be chosen 
BO provide a measure of those things that are of interest in 
Ewparticular system when it is used for a particular purpose. 

Although a general statement of the useful forms of the 
Cost functions for switched-parameter systems has not yet 
been made, there are some indications that the full gener- 
eilaty of (24) and (25) should be restricted depending upon 
MrcevariOus applications. If the control interval sb fone 
fimea length for a particular problem, the weighting factor, 
R, may as well be zero since the RT term in (24) will then 
only add a constant to Jin nea part ictm@er<ubinverval 


th 


eae i , say) the matrix M. would have non-zero elements 


for each state whose deviation from a specified trajectory 
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was to be penalized (i.e. increase the Cost function value). 
Likewise, Ny would have non-zero elements for those elements 
of the control vector, u(t), whose scaled magnitude was to 
contribute to the cost in the subinterval. A similar state- 
ment applies to He and the deviation of the values of the 


state variables at the endpoint of the subinterval from 


D 


Variation of a system variable is to be ignored, the cor- 


some specified goal, x (t.,4)- Of course, if the value or 


responding element in the appropriate weighting factor 
imeiwier x 1S Set to zero. 

The Cost function measure of system quality may be used 
in a function minimization scheme to optimize the variable 
me eters in the. system. im so doinge@yeeone Urears Jin as a 


P FWA wWwANM 
he ae ae Nae ee 


Hom C Fs 
function to obtain "optimum" values of the parameters. For 
example, one might optimize the switching instants (the 
ts for i= 1,N) for a particular switched-parameter system 
(er a Particular input, or one might try to minimize the 
total control interval necessary for the system to complete 
its response to another input (by Sel sata varying the t,, 
Suge initial condition due to an isolated state, perhaps), 
or attempt to optimize one of the system parameters. 

In general, the possible values of the independent vari- 
ables of the minimization will be limited (e.g. the Us may 
be constrained such that 0 < t, < . fOmme ct Genel). aid 


toe Minimization of Jin must use those procedures applicable 


to constrained systems. It may be noted, in passing, that 


Ded 





the control vector for the switched-parameter system may 
be chosen so as to Minimize @ fumerion aio! Se temo diieen Jims 
This represents the optimal control approach to determining 
mae control function for switched-parameter systems, it 
is an area that remains to be explored. 

As examples of their use, two Cost functions were con- 


structed for the system of (19) and (20) and Figure 2. The 


iomest Cost function was 


: aL 
x, (T,) = 116 0 0 0 x4 (T, + © 
Jn = oe be 0 il 0 "Sy ae IL ee 
x3(T,) 0 0 0 1 X2(T,) - 0 


Jn was treated as a function of the length of the first 


ote 3 -. 2&- w---+- 7 4+ A, sm = 
Sie wae ia ¥ ws oi va ° 


c “> <i - @ w- 4 = - 7 din G are : 4. J - eK - Gi) bo 
, OUGe Oma Sued oi .Oi uhe sy.cec etree 


(19) and (20), over the control interval AN 5 evens eval- 


Matton of (26) for 100 incremental values of t. spanning 


1 
the interval lo,2] , wilt a = LOe voll Ste) fPuncelcom ea eo 
the system produced the points for the curve plotted in 


Figure 8. 


fienires 9 and 10 Show the Cost function of (26) modified 


so that deviation from the desired final value during the 
entire second SUB Piive “creas counted as part of the cost. 


iiemmodified Cost function is, then, 


a 








20 


12 


Oot 


~ 
- 


Figure 8. 


fl 





0.4 0.8 ee 1.6 2.0 
TIME TO FIRST SWITCHING IN SECONDS 


Coop PuMmeulen Om cf tit, GrOveeciam ol Lened—raramnever 
System as a Function of the Time to the Farst 
Switching for a 10 Volt Step Input and Where Only 
the Square of the Deviation of the States from 

the Desired Final Value is Counted as Cost. 


54 








12 


3 
rae 
Oo 
x 
| 
Vv) 
O 
aS 
Figure 9. 


= 20 


15 


/ 


© 
\ ; / foi 
3 f : 
9) 
O 0.4 0.8 he 1.6 2-0 


TIME TO FIRST SWITCHING IN SECONDS 


The Cost Function of Figure 8 Modified to Count 
(ne Limominvesralsor the sduare of view Devialion 
fremmcne Desired Hinal State During tae Entire 
Secona COnumol "Subinverval 2s CoOst,etor Syovem 
lippuersvCDs OL Ver1eus Ape wudes:. 


D> 





© 
20 


QO 
Ty; 
= 
6 Vv) i 
S- e 
O a 
ce esi) ak 
= 
WY) 
O 
O) 4 l-. e ea. 
2 
2 
A 
O 1 2 3 4 5 
MIME TO.FIRSt sw CeCHING IN=sEeeCOnes 
Figure 10. Same as Figure 9 Except that Time to the ILO 
Swicehine Te Aivowed tomvary in the Interval 


[0,5]. 


56 








x, (T,) - 0 0 0 0 x, (T,) - 0 
Jin a ae 0 a 0 a aaa + 
X(T) - 0 0 0 ili X(T) ) - 0 
(27) 
qT. _ 

x, (s) - 0 0 0 0 x, (s) - 0 
X5(s)- XxX, 0 a 0 X5(S)- xX, ds 
X2(s) _ °| 0 0 0 X2(s) -~ 0 

ci 


The scalar x, is both the desired trajectory (a constant) 


and the desired endpoint value for the second subinterval. 
ites also the amplitude of the input step function. 
The purpose of (26) and the first term of (27) is the 


same; namely, to penalize (increase the value of the Cost 


-, Bel ey Ph Lo, ~ = he = + AMR tT 
i ee pe - 


Mrevid:n, Vie GSVWation of ti Aves 
Vee s at the endpoint of the control interval. The Cost 
mee 1 On, Jims is increased by the square of such deviation 
as does exist at that time. In both (26) and (27) the 
values of x5(T.) and X(T.) are to be penalized if they 
Gagm@ner from X, or zero respectively. The endpoint value 
of x, (t) does not affect the value of the Cost function 
Since the first element on the diagonal of the weighting 
fewor matrix, i, is zero. 

The second term of (27) serves to penalize the system 
gimethe trajectories of its states vary from some prescribed 


Heeyectory during (in this case) the second subinterval. 


Since the only nonzero element of M corresponds to the 


1 


deviation of x5 (t) from x only the integral of the square 


qa? 
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of that deviation will add to the Cost function. The result 
is that both overshoot and slow rise time of Xo (t ) in the 
second subinterval for the switched-parameter system being 
considered will be penalized. 

mf (27) were tempe applied to a shaft positioneeencrol 
system with constant parameters, and if Xo (t) and X(t) were 
the state variables corresponding to the shaft position angle 
and its time derivative respectively, then the application 
of the Cost function may be seen to give a measure of how 
emosely the system approaches the value of the input step 
fmec ion after some period of time which is considered as 
woewconvrol interval. The Cost function would also give 
a measure of how rapidly, without producing excessive over- 
Paseo Or Lie wena Var , vir Sysvem ScCus vo tnace Finel 
value. It may be seen from consideration of this example 
That the Cost function is ambiguous in that an overdamped 
system which rises slowly to its final value may have the 
same value of cost as does an underdamped system which over- 
shoots and then oscillates about the final value several 
Limes. 

The multiple curves in Figures 9 and 10 show how The 
Meer, function varies as the amplitude of the input step in 
varied. Figure 9 is analogous to Figure 8 in that t, is 
oes allowed to vary in the interval [0,2]. Figure 10 
Saews how rapidly the Cost function increases as vy is 


increased beyond 2 and provides a non-rigorous basis for 


believing that the minimums evidenced by the functions in 
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Figure 9 are truly global minimums. The system's states 
were all started from Zero initial conditions Giaeu ie colon. 
tation of these curves. The possibility of minimizing the 
Coste functions of Figures S and 9 with respect yo ne 
"Variable" ty is apparent. ; 

The determination of Cost function plots such as Figures 
8, 9, and 10 requires a large expenditure of computer time 
since repeated solutions of the state equations are neces- 
sary to evaluate the Cost function for different values of 
its independent variables. Each individual curve in the 
three figures required approximately three minutes of compu- 
poawen time on the IBM 360/67 at the Naval Postgraduate 
Semeol's W.R. Church Computer Center. These plots are for 


ey an isto gone oat eno 
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Wem Mave been required for more than ome independent 
WVemmable. Fortunately, function minimization usually 
requires fewer evaluations of the Cost function and, hence, 


less computation time. 
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V. SUMMARY AND CONCLUSIONS 


The preceding material has 1) developed a general 
description of linear, switched-parameter systems in terms 
of the mathematics of state variables, 2) presented the 
concepts of state variable grouping and of the "core" 
states of a switched-parameter system, 3) presented two 
procedures for determining the sequence of coefficient 
matrices in the state variable representation That describes 
the behavior of the switched parameters in the system and 
mim@osurated their use with simple examples, and 4) presented 
one possible Cost function which may be used to measure 
Unemer i ormance cf those @yctems ermaeannites if «mn two 
examples to demonstrate its application as a measure of 
Soeem Quality. That such a Cost function may provide a 
useful way to apt Sa ae the behavior of the system has been 
demonstrated in Appendix A. 

Several conclusions may be drawn from the work just 
Summarized and the examples which were presented in conjunc- 
tion with that work. Although the model that was developed 
was not applied to a real-world system, it is obvious from 
both the mathematics and the examples that physical systems 
like those presented in the Introduction may, indeed, be 
modeled as Switched-parameter systems. A conclusion that 
is possible as a by-product of the examples of switched- 


parameter systems that were used is that it is possible to 
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greatly improve the step response of a heavily damped second 
order system by using the states in the second order svsren 
as. "core" states in an unstable second or ¢hivd order 

system in order to "pre-charge" them. The penalty that one 
must pay for such improved response is a marked dependence 
of the response on the time that is taken to "pre-charge" 
the states. 

One may also conclude that it is possible to use Cost 
function minimization to. determine the optimum switching 
time which produces the "best" step function response in 
at least one variety of switched-parameter system. The 
optimum switching time so determined is dependent upon the 


form of the Cost function which was minimized. Finally, 


= Coe Pe aa me Sf Jo FH is 
lav Vv s11G MoStuat Ww Waipu 


tes evident tT Rac ctouMmsed sta 
solve the state equations and to determine the optimum 
Switching point were both sensitive to and able to cope with 
the sensitivity of the system response to the switching 
time. 

This first attempt to study linear, switched-parameter 
systems in some generality has raised many more questions 
than it has answered; indeed, much the greater part of a 
complete study of these systems remains to be done. Not. 
epm@eeot the possible paths that lead from this tentative 
beginning have been discovered. Some of those which have 


Meenmevand which may prove worthy of further study are pre- 


sented in the following section. 
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VI. RECOMMENDATIONS FOR FURTHER WORK 


The areas of further study that are recommended below 
constitute a study of what may be called the "design prob- 
lem." Some of the questions require that one attempt to 
answer the question, "What is the optimum system to use for 
this job?" That question has been answered heuristically 
ioe ne past and a definitive answer is not expected soon. 
Still, one may ask questions about the system architecture 
that do have some hope of being answered. 

iniewirirst problem rel@tange to the architecture of the 


systems that might be studied would be to define the con- 


NAC, OW tee G 4) Aaah) 


Stewie s On vite pirySicaly emctem hart arc reaniz 
eme-rrain order of behavior of the “core” states. For 
iiemeamee, one might ask what system architecture is required 
to produce third order behavior in a system with two "core" 
Seenve S . 

The kind of behavior that one might seek to establish 
in a switched-parameter system in order to satisfy the design 
goals should be related to what is possible, desirable and 
obtainable in the: way of system architecture. This. might 
be approached by seeking to determine what sort of output 
responses are possible when one switches between two par- 
ticular linear systems. If one were to form a table of 


systems before and after switching, and the particular 


response characteristics obtained, one might enter the system 
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used in Figure 3, noting that it 20es froneaee ede. der 
system with a pole in the right half of the s-plane To a 
damped second order system and that the result is an im- 
provement in the system step response for a suitable 
selection of switching time. The hope is that [by looking 
at such a table], one would be able to decide which system 
GO switch to in order to obtain a desired response. 

A “building block” approach may also be possible. From 
toes table just mentioned it may be possible to identify 
certain sequences of linear systems which have character- 
istics such that one might desire to use these subsequences 
as building blocks in constructing various systems. Of 
course, this raises the question of determining what sorts 
OF POULLGITE OLOCKS TM Ent De Tequirec, or cestred, T 
eeeoyscvem and how to put them together to obtain the desired 
Pesulls . 

The work of Flugge-Lotz and Taylor [6] that was mentioned 
in the introduction indicates that a quasi-adaptive scheme 
can yield useful results. A generalization of those results 
in terms of switched-parameter systems would be desirable. 

For switched-parameter systems generally, one might 
inquire what the switching criteria should be. As an example 
one may not want to switch to a particular system while a 
capacitor has a non-zero charge. One might switch instead 
(with the switching based on the values of some other states) 
to an intermediate system which would discharge the capaci- 
wer bDetoOrewmaking Whe transition to the articular system 


firsts. mene LOnec . 
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As a grand goal of the work with Syst emearen wee unre 
one might seek a procedure for the systematic design of a 
switched-parameter system to meet some specified require- 
ments. 

Optimization problems represent another area of study 
that could produce significant results. One usually encoun- 
wemgsesuch problems in connection with optimizing an input 
to an already fixed system; however, with a system which is 
moe fixed (e.g. switched-parameter) it may be reasonable 
to ask what sequence of systems gives an optimum response 
iemea fixed inpuG. This has the nature of an inside-out 
optimal control problem where one seeks to determine the 
system that makes a given control optimal. 


be alle ohn) --+ cae +L 
vu 


me possivLliLey OF eo RSC ee ee a). 
certain of the states in the system was shown in Section III 
to allow a large improvement in the response of a heavily 
damped second order system, and has been mentioned as a 
Meeetble candidate for an optimization technique. This 

idea should be explored further, as should the theoretical 
improvement in system response that may be attained. Being 
oeeee CO vary the initial condition of one of a system's 
Staves would seem to permit one to use the transient response 
ig@ecne System caused by jGhat. initial condition to alter the 
transient response ofthe other states of the system in a 
@om@erolled fashion. 


ihe land Of COSE IUnmeelon thae One Usecomwidl inttuenice 


the results of the optimization of a system. The 
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correspondence between changes in the Cost function and 
changes in the physical response of the system is very 
important. As was seen in the various examples, a smooth, 
slowly varying Cost function may correspond to a physical 
mowwense which 18 extremely Sensitive to The panrueiiar 


parameter variation, as was the case with Tt The estab- 


1: 
lishment of guidelines for the assignment of weighting 
factor values, thus, seems useful. It may be worthwhile 
To investigate the type of Cost function necessary to allow 
a system parameter to be optimized with the result that 
some of the system's states may be constrained to stay 
within certain limits. This would be useful since it would 
allow the designer to specify physical limits on the system's 
states where such limits resulted rrom saturation iimits on 
the physical device or from some desired performance criteria. 

items iOuld oe possible to formulate an Optimal control 
problem for switched parameter systems. This should be 
@emwermined and any limitations on the system imposed by the 
requirements for the existance and uniqueness of the optimal 
control systems with constant parameters and with the "opti- 
mal" control system mentioned previously in which the 
femerol 15 obtained by switching parameters. 

ire descriving the Switching curves of the present relay 
femerol SYSteMmS, it 1S common practice to speak of the 
Switching boundary as a hyper-surface in state space. Fur- 


ther work might be done to determine the effects (beneficial 


or deleterious) of being able to switch on a particular 
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"line" in that hyper-surface, where such flexibility would 
be allowed by the independent contro! of Chewinitt saeco waa 
tion of a particular state which was isolated from the 
"core" states at some point. Such considerations could 


also be the basis for the design of certain small segments 


Soma larger system. 
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APPENDIX A 


APPLICATION OF COST FUNCTION MINIMIZATION 


The purpose of this appendix is to demonstrate the 
utility of a Cost function minimization procedure as a 
means of optimizing the step function responses of some 
Simple switched-parameter systems. As the first example 
eeearcOSt function to be minimized, consider Figure iI2. 

The set of functions shown there resulted from the appli- 
cation of five different amplitude step functions to the 
Switched-parameter system of Figure 11, which is the system 
of Figure 2 reproduced here for convenience. The system 
ferme cen 1ruiln Zerc inivial ccnaitions for ead eefates in 


each case. 


u(t) 





me; 'T' | 


Breure J1. Replica of Figure 2. Block Diagram of a Third 
Order, Switched-Parameter System. 
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Bunciaton Of apne Tame toe the First switching for 
Various Amplitude Step Inputs; where the Square 
Of Bothwthne, Deviation ot the sraves trom the 


Desired Endpoint Values and the Deviation of 
the States from the Desired Final State During 


the Second Subinterval are Counted as Cost. 
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From Figure l2 it is easily seem thav theres stumer lene 


Boch have a2 Minimum fOr a time vould tsi seu elie ere of 


1? 
epout 1.35 seconds. For low amplitude steps (e768. 5 volte 
or 7 volt) the minimum is quite broad. It becomes more 
sharp as the step amplitude increases to 20 volts. Because 
of the broad minimum, Qne would expect the time eo ae 
@eecne system not to vary drastically for a small change in 


18 A comparison of Figures (13(a) and (b) (which are 


ples 
mepamecas Of Figures 3(a) and (b)) shows that this is the 
case. The small increase in the Cost function which is 
barely discernable on the 10 volt input step curve of 


feewiee 12 as t. increases to 1.45 seconds, i8 seen to ve 


1 


reflected in an overshoot-and-decay response in Figure 13(a). 
| S the Yaect thet Lhe Various 2s 
functions of Figure 12 all have about the same minimum. 
Hemeen lor this system, if the switching time is taken as 
Mmeooeopeconds, the system response should be near optimum 
/memea wide variation of step inputs. The system behaves 

eas a linear system, the character of its output response 
being unaffected by the amplitude of the input. ‘This linear 
response is to be expected since the system is linear in 
both subintervals and the switching instant, Cy» is inde- 
pendent of the values of the states, being fixed in value. 
Mitemis brought out vividly in Figures 14, 15, and 16, which 
show the output response of this system for input steps of 


peer, and 15 volts respectively. 
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These figures should be compared with Figures 13(a) and 
(ob). Note that there is 4a slight ehange agneeher ve uemom 


merfrom 1.349 in Figure 13(b) to. 173475 torebwetreos) Ge 


i) 
and 16. 

ine second, example of the appdiear om Oto e fea epake nn 
minimization will be for a second order system with a 
Switched damping factor. The block diagram of the system 


before and after switching is shown in Figure 17. The 


state equations for this system are 


: Oe 1 0 
(rie Faron: (ular eK 
- a Onl 1 
te[0,t,], 
(a) te[0,t,] 





(b) te[t,,10] 


peeure lf. Block Diagram of a Second Order, Switched-— 
Parameter System. (a) Before Switching. 
(bo) After Switching. 
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and 


mt) = x(t) + [ues T sen 
hg Tae il 


te[t,,10], 


where the length of the control interval is taken as 10 
eee@onds. 

A Cost function for this system was constructed analgous 
to the one for the third order system just discussed, the 
PUimerihor being dependent on the final values of the states 
mumeene end of the control interval and dependent upon the 
deviation of the state x, (t) from the desired final value 


during the second subinterval. Thus, the Cost function was 


7 Px, (20)-x] [2 O}[ x4 (10)-x,] 
FL x,(20)-0 | | o 1][-x5(10)-0 | 


hg 


(28) 


x, (s)-x, 7 al 0 X,(S)-X, 7 
X5(s) -0 0 0 CD ae 
ei 

The first term of (28) penalizes the deviation of the 
Staves from their desired values at the end of the control 
iiewermval by increasing the Cost function by the sum of the. 
Squares of the deviations of both SEaCee from the desired 
end point values. Deviation of x, (t) from the input step 
amplitude (the desired final state for x4) represented by 


Xa. and X5(t) from zero final amplitude are both penalized 


in this way. The second term penalizes any deviation of 


15 











x, (t) from the desired trajectory (in this case, a constant 
value, Xa which is also the desired final state) during the 
second subinterval, by adding the integral of the square 
oemene deviation to the cost during the Subiyjverva boas. 
penalizing the deviation from the desired trajectory in this 
way it is hoped that the optimization procedure can be forced 
memechnoose a switching time, Uy > which will give the system 
weeaplad rise time and little overshoot. This Cost function 
iemolertcted as a function Olt the time to the first svieenmmnea. 
Cys in Figure 18. Zero initial conditions were imposed 
upon x, (t) and Xo (t) in depermiI nine thts eer funeume ne 

In order to test the usefulness of the function minimi- 
Bagmon approach, the Cost function, Jn> was minimized as a 
Pomezilon olf ty using a Fortran iv submoucrlWe with amp vemencs 
the Tangent Search method for the minimization of a function 
with constraints. See [9] and [10] for further details of 
the program and the method. 

[he result of the optimization is shown as Figure 19, 
which shows the "optimum" time response of this second 
Cee. Switched-parameter system as the functions x4 and Se 
These response curves should oe Compared TOrcniosce ©1 the 
constant-parameter second order system that are shown as 
xy and Xo which are taken from the response curves of 
Figure 4. The improvement is evident, as is the price that 
must be paid for "snappy" response; namely, that the time 


response of the velocity Cre) becomes impulse-like, indi- 


catTing that the system of x4 and Xo is undergoing some 
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relatively violent accelerations compared to the soft res- 
ponse of the constant-parameter, linear system. 

The final example demonstrates the improvement that can 
be made in the step response of the third order, switched- 
parameter system of Figures 5 and 6 (shown below as Figure 
20). The time responses of this system to a 10 volt step 
input, for two arbitrary switching times are shown in 
Figure 7 which is reproduced below as Figure 21. 

Two Cost functions were constructed for this system. 
Msetor the previously mentioned Cost functions, these two 
were plotted assuming that the system's states started at 
Zero initial conditions. One was dependent on the values 


Seeeeme states at the endpoint of the control interval only, 
Gupmeeeiceese (26) and PZstume © for thossrstem cf Figur: 22; 
and the second Cost function was analgous to (27) and 
Figures 9 and 10. The Cost function which depends only on 
Eiiewdeviations of the endpoint values,of the states from 
meomemeaesitred final value is shown in Figure 22 for a system 
femea 10 volt step input. It is seen to be not markedly 
different from Figure 12 for the 10 volt step curve except 
miaeeene Minimum value is displaced slightly to a larger 


value of t,. The Tinimumets seem compe forur siliehcly less 


il 
than 1.6 for an essentially zero value of the Cost function. 
Pilphoueh it is not evident from Figure 22, the Cost function 
Pouwinues its rapid rise as ty is increased beyond by = 220 


in the same manner as that shown in Figure 10, implying that 


the minimum at t, = i omes =a velop ad omit meine 


(ce 








(a) 





te[0,t,] 


(b) u(t) 





Figure 20. Replicas of Figures 5 and 6. Block Diagrams 
Cla eeedeO rice We Vehed— ra rwame lel aoy Sei, 
Caer lore rosso wirecning. Go) Aiter.) sowie champ. 
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This Cost function was minimized using the same method 
as was used for the second order system just considered with 
the addition that three different values of step input 
amplitude were used. The "optimum" switching time was found 
to be the same for all of the inputs, t, = lL. 562 seconds. 
The time responses of the output states, X(t), for the 
optimum response to the various step inputs are shown in 
Figure 23. Figure 24 shows a similar plot of the time 
responses of X(t), the derivatives of the output states. 
ict should be noted from Figure 23 that the character of 
Une Optimum response is independent of the input step 
amplitude, the various responses differing only in magni- 
Medammetnis result is not unexpected. It is a reflection 
eee J oeCn' s lineaMiuy in each subinterval emda the 
independence of ty of the values of the states during the 
montwel interval, just as the third order system considered 
at the beginning of this appendix showed linear behavior 
Metmine Same reasons. 

in order to obtain some feeling for the influence the 
meomess Of Cost function has on the “optimum” response 
©bGained, or, what is the same thing, on the value of by 
preeauiced by the minimization of the Cost function, emou mice 
West function was constructed for the same system. This 
new function was the analog of (27) and Figures 9 and 10. 
Miiemmew COSt function is plotted as a function of the time 


Mommene first switching in Figure 25. There, it may be 


noted that the "optimum" value of ty is not the same as 
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it was for Figure 22; but 1s increased mm vale se tinct aco 
ebout 2.6 seconds. 

Again, a function minimization was done on the Cost 
function. The "optimum" switching time so obtained was, 


1 


we. O00 seconds. This optimization was also acconplmened 
for the same three amplitudes of input step as were facia 
in Figures 23 and 24. The "optimum" trajectories, X5(t) 
and X(t), SremolOnced Line ie cme omanicenc ce 
The main difference to be noted between the responses 
peeeeweures 23 and 26 is that the first “optimum” trajectory 
shows a slight undershoot, and the second a slight over- 
oeoee,., Both optimum systems show a marked improvement over 
that shown by the same system in Figure 21. As was the case 
foeeenes trrs. Cost Tunetyvon considered Lor tise SVStSM. Bic 
optimum switching time is seen to be independent of the 
input step function amplitude, and for the same reasons. 
Cosy function Minimization is, thus, seen to providewa 
ieeamui means of “optimizing” this particular kind of 
switched-parameter system. The utility of this kind On 
procedure for improving the responses of other kinds of 


Switched parameter systems remains to be demonstrated. 
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